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Abstract 

o ■ 

^sj ' Let the Bessel number of the second kind B{n,k) be the number of set partitions of \n\ into k 

blocks of size one or two, and let the Bessel number of the first kind 6(n, k) be the coefficient of x^~^ 
in —yn-i[—x) , where yn{x) is the nth Bessel polynomial. In this paper, we show that Bessel numbers 
Q_^ ' satisfy two properties of Stirling numbers: The two kinds of Bessel numbers are related by inverse 

T\ . formulas, and both Bessel numbers of the first kind and the second kind form log-concave sequences. 

By constructing sign-reversing involutions, we prove the inverse formulas. We review Krattenthaler's 
injection for the log-concavity of Bessel numbers of the second kind, and give a new explicit injection for 



^ 



C^ ■ the log-concavity of signless Bessel numbers of the first kind. 

1 Introduction 

^ i 1.1 Main Results 

C*^ ' For nonncgative integer n and k , let B{n, k) be the set of partitions of [n] :={!,..., n) into k blocks of size 

one or two. Let B{n, k) denote the number of elements in B(n, k) . The Bessel polynomials 012 El E] are 
(*-^ ' the polynomial solutions y„ {x) of the second-order differential equations 

O; x'y'^ + {2x + 2)y'^^n{n+l)y,,, 2/„(0) = l. 

"^ . Let a{n,k) denote the coefficient of x"^'' in y„-i(a;). In fact, a{n,k) equals the cardinality of the set 

<^ ; B{2n -k-l,n-l). Set b{n, k) := (-1)""'= a(n, k) . 

In this paper we deduce following two identities. 



X 



J2B{n,k)b{k,l) = SnA and J^ b(n,k) B{k,l) = Sn,i ■ (1) 

fc=0 k=0 

These identities are reminiscent for the inverse formulas of Stirling numbers [S] Prop 1.4.1]. By constructing 
sign-reversing involutions, we give combinatorial proofs for both identities. 
Second, we show the log-concavity of { a{n, k) }i<fe<n , i.e., that 

a{n, k — 1) ■ a(n, A; + 1) < a{n, k) , 

which is an analogue of the log-concavity of the signless Stirling numbers of the first kind. We prove this 
result by constructing an explicit injection from the set B{2n — k,n — 1) x B{2n — fc — 2, n — 1) to the set 
B{2n — k — l,n— 1) x B{2n — /c — 1, n — 1) , which provides a combinatorial proof of the log-concavity of the 
sequence. 
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1.2 Notations and Backgrounds 

We review some standard terminology, which wiU be used throughout the paper. 

Stirling numbers 

For a nonnegative integer n , let [n] denote the set {1, 2, . . . , n} , and let [0] = . A partition of [n] is a 
collection of nonempty pairwise disjoint subsets whose union is [n] . These disjoint subsets of [n] are called 
blocks of the set partition. The number of partitions of [n] with k blocks i?i , . . . , Bk is called the Stirling 
number of the second kind S{n, k) . By convention, we put 5(0, 0) = 1 . We denote by Sn the set of all 
permutations of [n] . Let c(n, k) be the number of permutations in Sn that have k cycles. Define s{n, k) to 
be (—1)""'^ c{n^ k) . This number s{n, k) is called the Stirling number of the first kind. It is well-known that 
Stirling numbers are related to one another by the following inverse formulas 9, Prop 1.4.1]: 

n n 

'^ S{n,k)s{k,l) = dn,i and ^^ s{n,k) S{k,l) ^ Sn.i ■ (2) 

fc=0 k=0 

Moreover, both Stirling numbers satisfy the log-concavity property, i.e., that 

cin, k - 1) ■ cin, k + 1) < c{n, kf and S{n,k - 1) ■ S{n,k + 1) < S{n,kf . (3) 

There are many proofs for (01 and Q with both algebraic and combinatorial ways. See pi 1101 ITT] . 

Bessel numbers 

Recall that the Bessel polynomials are the polynomial solutions y,i{x) of the second-order differential equa- 
tions 

x^y'; + {2x + 2)y'^ = n{n+l)yn (4) 

satisfying the initial condition y„(0) = 1 . Bochner pP seems first to have realized that these polynomials are 
closely related to the Bessel functions. Krall and Frink "W have considered the system of Bessel polynomials 
in connection with certain solutions of the wave equations. Moreover, Burchnall 2, has developed more 
properties of the Bessel polynomials in detail. Recently, Choi and Smith 3^ have investigated coefficients of 
Bessel polynomials from a combinatorial point of view. 

From the differential equation Q one can easily derive the formula: 

Let a{n, k) denote the coefficient of x^^^ in yn-i{x) . Set 6(n, k) := (—1)"^'^ a{n, k) , i.e., 

b{n,k):=l^ ' 2^-^{n-k)\{k-l)\' " ^ - '« - " ' (5) 

[ , if 1 < n < fc . 

We call the number b{n, k) a Bessel number of the first kind, and a{n, k) a signless Bessel number of the first 
kind. By convention, we put a(0, k) = 6(0, k) = (5o,fc • 

For nonnegative integers n and k, define B(n, k) to be the set of partitions of [n] into k nonempty blocks 
of each size one or two. A block B is called a singleton ii \B\ = 1 , and a pair if \B\ — 2 . Note that if 
TT £ B{n, k) , then n has exactly 2k — n singletons and n — k pairs. Set B(n, k) :— \B{n, k)\ . To choose n — k 
pairs, we should choose 2n — 2k elements from [n] and pairing the chosen 2n — 2k elements, so the number 
_B(n, k) is given by 

!n\ 
2-Hn-k)\{2k-n)\' ^^ T^/Sl < fc < n , ^^^ 

, otherwise. 

We call the number B{n, k) a Bessel number of the second kind. 

Remark 1. It is easily checked that a{n,k) — B{2n — fc — l,n — 1) . So we can consider a{n,k) as the 
cardinahty of the set B{2n — k — l,n — 1) . 



Matching number 

Let G be a loopless graph with n vertices. A k-matching in a graph G is a set of k edges, no two of which 
have a vertex in common. A matching a saturates a vertex a; , or x is said to be saturated under a , if an edge 
in a is incident with x ; otherwise, x is unsaturated under a . The complete graph Kn with n vertices is a 
simple graph with each pair of whose distinct vertices is adjacent. Let A4{n, k) denote the set of fc-matchings 
in complete graph K^ and m{n, k) be the cardinality of A^(n, k) . We call m{n^ k) a matching number. From 
now on, assume that all vertices of K^ are labeled with the integers in [n] . Under this assumption, a partition 
of [n] into k nonempty blocks, each of size at most two, can be identified with an (n — fc)-matching in Kn 
canonically. For example, the partition tt — { {1}, {2, 6}, {3}, {4, 5} } in S(6, 4) , corresponds to the matching 
a = { {2, 6}, {4, 5} } in M{6, 2) . Hence B{n, k) = m{n, n — k) . Whether {a, 6} is an edge or a pair, we will 
write it such that a < b . 

It is well known (see ^ElCj) that the matching numbers { m(n, k) }o<k<n form a log-concave sequence, 
i.e., that 

m{n, k — 1) ■ m{n, fc + 1) < m{n, k) , (7) 

which implies that the Bessel numbers of the second kind B(n,k) is also log-concave. In particular, Krat- 
tenthaler ^ has given a neat combinatorial proof of 0. We review Krattenthaler's proof in Section EITl 

2 Inverse formulas for the Bessel numbers 

In this section we deduce the inverse formulas for the Bessel numbers Q by the generating function technic. 
For indeterminates x and t, let 



n=0 



t" 



l + t + 



f \^ 



2! 



Note that { fn{x) } is a sequence of polynomials of binomial type, i.e.. 



A:=0 



fnix + y)=Y^ ( ',_" ) fk{x) fn-k{y) ■ 



For details on polynomials of binomial type, see |H1 pp. 167-213]. 
Lemma 2.1. For any nonnegative integer n, we have 

n 

fn{x) =^ B{n,k) [x]k, 

k=0 

where [x]q — 1 and [x]k ~ x{x — 1) • • • (a: — fc + 1) , for fc > 1 . 
Proof. We have 



fn{x) 



f' 



t' 



^l(l+t+^ 



E 

fc>0 



x\ I t'^ \k 
fc (^+2! 






fe>0 



jt + kr 

fc! 



[x]k 



= Y,B{n,k)-[x]k. 



fc>0 



The last equation holds due to the exponential formula ^2 p. 81] . 



D 



Lemma 2.2. For any nonnegative integer n , we have 

n 

Nn =X! Kn,k) Ik{x). 

fe=0 

Proof. From the formula © and Lemma [2. 11 we have 
Multiplying both sides of this equation by T^ /n\ yields 



fc=0 



fc \ 2'= 
n — k I k\ 



^ L_ ,J — Mfe 



Wilf 11, p. 169] shows, using the Lagrange Inversfon Formula, that if two sequences {a„} and {&«} are 
related by 

then we have the inversion 



eC";- ;')'-''""■ "^ 



k 

So we have 

Multiplying both sides of this equation by 2^" (n — 1)! yields that 

From the formula (O, we have 

[x\n = ^ b{n,k) fk{x) , 

which completes the proof. D 

Note that both [x]n and fn{x) are monic polynomials of degree n. By Lemma [2.11 and U^ we have the 
following inverse formulas for the Bessel numbers. 

Theorem 2.3. Let b{n, k) he the Bessel number of the first kind and B{n, k) be the Bessel number of the 
second kind. Then two Bessel numbers are related to as follows: 



Y, B{n,k)h{k,l) = ^„,z, (8) 

fc=0 
n 

Y,Kn,k)B{k,l) = S^j. (9) 

fc=0 

In Sectional we prove these two inverse formulas © and © combinatorially. 





Figure 1: The involution /i 



3 Two involutions for the inverse formulas 

Reeall that a partition of [n] into k blocks of size one or two can be identified with an (n — /c)-niatching in 
Kn- So we can describe the Bessel numbers in terms of the matching numbers m{n, k) as follows: 

B{n, k) = m(n, n - k) and bin, k) = (-1)""'= a(n, k) = (-l)""*^ m{2n -k-l,n-k). 

3.1 An involution for the first inverse formula 
Theorem 3.1. For all nonnegative integers n and I , we have 



^B{n,k)b{k,l) 



fe=0 

Proof. Consider the complete graph K2n-i-i- For arbitrary k , let U be the set of ordered pairs (a, (3) such 
that 

1. a is an {n — fc)-matching in K2n-i-i whose saturated vertices are in [n] ; and 

2. /3 is a (fc — /)-matching in K2n~i-i such that a and (3 saturate no vertex in common. 

Define the sign of an element (a,/3) G C/ by sgn(a,/3) = (— l)!''' , where |/3| denotes the number of edges in 
(3. Then we have 

Y^ sgn(a, f3)^Yl (-1)''"' »"(". n-k) m{2k -l-l,k-l) = Y^ B{n, k) b{k, I) . 

{a,l3}eU k k 

To prove the theorem, it suffices to find a sign-reversing involution /i on U having exactly one fixed point 
with sign +1 . 

Consider the union of a and /3 . Define a linear order on a U /3 as follows: For edges {a, 6} and {c, d} in 
aU (3 , {a, b} < {c, d} if and only if 6 < d . Note that this order is the same as the colex order (see ^1 p. 8]). 
With respect to this order, take the smallest edge e of a U /3 and move e from a to /3 , if e is in a ; otherwise, 
move e from /3 to a . This defines the needed involution /i on [/ . 

Since there are n — I edges in a U /3 , for the smallest edge {a, 6} in a U /? , the number of vertices greater 
than b in [2n — Z — 1] is at least n — I — 1 . This implies that 6 < (2n — I — 1) — {n — I — 1) = n . So /i is 
well-defined. Obviously, Ji changes the number of elements of (3 by one, i.e., /i is sign-reversing and double 
applications of /i on (a, /3) will clearly restore (a, /?) . This construction cannot be accomplished, if both a 
and P are empty sets, i.e., n = k = I . Thus we obtain the desired result. D 

Example 1. Given n = 7, /c = 5 , / = 2 , let a = { {2, 3}, {4, 7} } and /3 = { {1, 10}, {5, 11}, {8, 9} }. 
Then a and (3 are matchings in Ku satisfying (a,/3) G U . Take the smallest edge e = {2, 3} in a U /?. 
Since e is an element of a , the image of (a,/3) under Ji is {a',/3') , where a' = a \ {e} — { {4, 7} } and 
/?' = /3 U {e} = { {1, 10}, {2, 3}, {5, 11}, {8, 9} } . It is easy to check that /i(a', (3') = (a, P) . (See FigureH]) 




< — - — > 




Figure 2: The involution I2 ■ 
3.2 An involution for the second inverse formula 

Theorem 3.2. For all nonnegative integers n and I , we have 

n 

Y^b{n,k)B{k,l)^5r.,i. 

Proof. Consider the complete graph K2n-k ■ Let Vk be the set of ordered pairs {a, (3) such that 

1. a is an [n — fc)-matching in K2n-k in which the vertex 2n — fc is unsaturated under a ; and 

2. /3 is a (fc — ^)-niatching in K2n-k such that a and /3 saturate no vertex in common. 

Set V — [jVk ■ Define the sign of an element {a, (3) &V hy sgn(a, 0) = (— l)'"' . Then we get 

Y^ sgn(a, 13) = Y^ (-1)"-'^ m{2n -k-l,n^k) m{k, k-l) = Y^ bin, k) Bik, I) . 

{a,p)£V k k 

To prove the theorem, it suffices to find a sign-reversing involution I2 on V , which has exactly one fixed 
point with sign -|-f . 

Suppose that (a, f3) €Vk ■ Now consider the largest edge e of a U /? according to the colex order in a U /3 . 
Define the map I2 as follows: 

• 1 f e e a , then set 12(0, (3) — {a\ {e} , /3 U {e} ) . Since e £ a , neither a nor f3 saturates the vertex 
2n — k . Moreover, the matching a \ {e} cannot saturate the vertex 2n — fc — f. So(a\ {e} , P U {e} ) 
belongs to 14+ 1 . 

• If e e /3 , then set hia, f3) — {aU {e} , f3 \ {e} ) . Regard ( a U {e} , 13 \ {e} ) as a pair of matchings in 
K2n-k+i ■ Clearly, the matching aUJe} cannot saturate the vertex 2n — fc + 1 . Thus (aUJe} , /3\{e} ) 
belongs to Vk^i ■ 

Evidently {I2 o l2){a,(3) = [a, (3) and I2 is a sign-reversing map, whenever a U /3 7^ . The unique fixed 
point of I2 is (a, /?) , where a and /3 are both empty, i.e., n = k = I . This completes the proof. D 

Example 2. Given n = 10 , fc = 8 , ^ = 5 , let a = { {2, 3}, {4, 11} } and /3 = { {1, 7}, {5, 10}, {8,9} }. Then 
a and (3 are matchings in K12 satisfying [a, (3) G Vs . Take the largest edge e = {4, 11} in a U /3 . Since 
e is an element of a , the image of (a, (3) under I2 is (a', f3') G Vg , where a' — a\ {e} — { {2, 3} } and 
/3' = /3 U {e} = { {1, 7}, {4, 11}, {5, 10}, {8, 9} } . It is easy to check that hia' , (3') = (a, /?) . (See Figure^) 



4 Log-concavity of the Bessel numbers 

In this section, we show that both Bessel numbers of the second kind B{n, fc) and signless Bessel numbers 
of the first kind a(n, fc) form log-concave sequences. In particular, we will check each result by constructing 
an explicit injection. 
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Figure 3: The injection Ik- 



4.1 Log-concavity of the Bessel numbers of the second kind 

Theorem 4.1. The Bessel numbers of the second kind { B{n, k) }fe>o form a log-concave sequence: 

B{n, fc - 1) • B{n, fc + 1) < B{n, kf for all k > 1 . 

Proof. In fact, this theorem is a special case of the log-concavity of the sequence of matching numbers W,^. 
For reader's convenience, we sketch here Krathenttaler's idea ^. Since B{n, k) — m{n, n — fc) , it suffices to 
show that 

m{n, n — k + 1) ■ m{n, n — k — 1) < m{n, n — k) . 

We will construct an injection Ik from the set of all pairs (ai, 0:2) of an (n — fc + l)-matching ai in Kn and 
an (n — fc — l)-matching 02 in Kn into the set of all pairs (/3i,/92) of an (n — fc)-matching /3i in Kn and an 
(n — fc)-matching P2 in Kn ■ Let an (ri — fc + l)-matching ai in Kn and an (n — fc — l)-matching a2 in Kn 
be given. Color the edges of ai and a2 in blue and red, respectively. 

Let G = G[ai Ua2] denote the subgraph of Kn induced by all edges of ai and a2 . Consider the connected 
components of G . Since ai and a2 are matchings, each connected component of G is either a cycle or a 
path. Note that all edges in the components of G are colored alternately. So all cycles and paths of even 
length have an equal number of blue and red edges. Meanwhile, in case a path of odd length, the number 
of blue edges in the path differs from the number of red edges by one. We call a path of odd length a blue 
path, if it has more blue edges than red edges, and a red path, otherwise. 

Let R be the set of all red paths of G and let r denote the cardinality of R . Since G has two more blue 
edges than red, the number of all blue paths of G should be r + 2 . Consider the set C of all blue and red 
paths of G . Label the paths of C from 1 to 2r + 2 in the increasing order of the largest vertices in the paths. 

Let A — {fli, . . . , Or} C [2r + 2] , where oi < • • • < a^ . Define (j){A) — AU {aj + 1} , where j is the largest 
i for which a^ — 2i is minimal (assume oq = 0). Then the map is an injection from all r-element subsets 
of [2r + 2] into (r + l)-element subsets of [2r + 2] . (See TIT, pp. 33-39 and p. 54].) 

Thus given the r-element subset S of [2r + 2] corresponding to _R, we have an {r + l)-element subset 
(j>{S) = 5 U {t}. Note that the element t corresponds to a blue path, say P. If we exchange the colors in P, 
then P becomes a red path, say P' . Now we consider the graph G' ^ {G \ P) L) P' . Obviously, G" has n — k 
blue edges and n — k red edges. 

Finally set /3i be the blue edges in G" and (32 be the red edges in G". Let //<-(ai,a2) — (/3i,/32) ■ It is 
easy to check that this defines the desired injection. D 

Example 3. Let ai = { {1, 2}, {3,4} {6, 11} {7, 12}, {13, 14}, {16, 17}, {19, 20}, {21, 22}, {23, 24}} and az = 
{ {2, 3}, {6, 7}, {8, 9}, {11, 12}, {14, 15}, {16, 21}, {19, 24} } . Then ai and 02 are matchings in K25 ■ Color 
the edges of ai in blue (— ) and the edges of a2 in red (••■). There are three blue paths 1 — 2 ■• -3 — 4, 
17 - 16 • ■ • 21 - 22 and 20 -- 19 • • • 24 - 23 and one red path 8 • • ■ 9 , which are labeled by 1, 3, 4 and 2 , 
respectively. So we have r = 1 and 5 == {2} . Since 0({2}) = {2, 3} , the blue path P = 17 - 16 • • • 21 - 22 
is changed to the red path P' = 17 • • • 16 — 21 • • • 22 . From the graph G' ^ {G \ P) U P' , we can extract 
two 8-matchings in K25 : Pi = { {1, 2}, {3,4} {6, 11} {7, 12}, {13, 14}, {16, 21}, {19, 20}, {23, 24} } and /Js = 
{ {2, 3}, {6, 7}, {8, 9}, {11, 12}, {14, 15}, {16, 17}, {19, 24}, {21, 22} } . (See FigureEl) 



4.2 Log-concavity of the signless Bessel numbers of the first kind 

Theorem 4.2. The sequence of signless Bessel numbers of the first kind { a{n, k) }fc>o is log-concave: 

a{n, k — 1) ■ a{n, fc + 1) < a{n, k) for all k > 1 . 

Proof. Since a{n,k) = B(2n — fc — l,n — 1) ~ m{2n — k ^ l,n — k) , this theorem is equivalent to the 
log-concavity of the sequence m{2n — fc — 1, n — fc) , for k = 0, . . . ,n . So it suffices to find an injection from 
the set of all pairs (^1,^2) of an (n — fc + l)-matching Ai in K2n-k and an (n — fc — l)-matching A2 in 
K2n-k-2 into the set of all pairs (i?i, B2) of an (n — fc)-matching Bi in K2n-k-i and an (n — fc)-matching 
B2 in K2n-k-i ■ 

Now we construct a new injection Is motivated by Krattenthaler's injection. Let an (n — fc -I- l)-matching 
Ai in K2n-k and an (n — fc — l)-matching A2 in K2n-k-2 be given. Assign the colors blue to the edges of 
Ai and red to the edges of A2 , respectively. Then there are two possible cases: 

(a) The vertex 2?! — fc is unsaturated under Ai . 

(b) The vertex 2n — fc is saturated under Ai . 

In case (a), if we disregard the vertex 2n — fc in K2n-k and regard K2n-k-2 as the subgraph of K2n-k-i 
induced by the vertices in [2n — fc — 2] , then Krattenthaler's injection Ik yields a blue (n — fc)-matching 
Bi in K2n-k-i and a red {n — fc)-matching B2 in K2n-k-i ■ If the vertex 2n — fc — 1 is saturated under 
Ai , then the blue path containing the vertex 2n — fc — 1 should be labeled with the largest number 2r + 2 , 
where r is the number of red paths in G[Ai U A2] . By the property of the map (f> ^J Chap. 2, ex. 25], the 
largest number 2r + 2 is not contained in {(f>{S) \ S) for any r-subset S of [2r + 2] . So a blue path containing 
the vertex 2n — fc — 1 is not selected to change into a red path. Moreover, the vertex 2n — fc — 1 cannot be 
saturated under A2 . Therefore B2 cannot saturate the vertex 2n — fc — 1 . 

In case (b), we introduce a new mapping as follows. Consider the sets X and Y of all unsaturated vertices 
under Ai and A2 , respectively. Let x be the vertex which is adjacent to the vertex 2n — k . Choose the 
element y in Y with the same relative order as a; in X U {x} . Since the cardinality of X U {x} is fc — 1 and the 
cardinality of Y is fc, the selection y is well-defined. Cut the blue edge Cb = {x, 2n — fc} and join two vertices 
y and 2n — fc — 1 by the red edge e^ . Then Ai \ {ei,} becomes a blue {n — fc)-matching Bi in K2n~k-i and 
A2 U {cr} becomes a red [n — fc)-matching B2 in K2n-k-i ■ Define In{Ai, A2) = (i?i, B2) ■ In this case, the 
vertex 2n — fc — 1 is saturated under B2 ■ 

Let Is be the map Ik 'J In, i-e., 

J /a'(Ai, ^2) , if Ai does not saturate 2n — k , 
Is[Ai, A2) = s 

\In{Ai,A2), if Ai saturates 2n — fc . 

Since both Ik and Ijq are injective and the saturating condition of the vertex 2n — fc — 1 makes that the 
image of Ik and the image of /jv are disjoint, the map Is is a desired injection. D 

Example 4. Let Ai = { {1, 2}, {3, 4}, {5, 10}, {6, 11}, {7, 12}, {13, 14}, {18, 19}, {20, 25}, {23, 24} } e X(25, 9) 
and A2 = { {2, 3}, {6, 7}, {8, 9}, {11, 12}, {14, 15}, {17, 18}, {19, 20} } e M{2i, 7) . Color the edges of Ai in 
blue (— ) and the edges of ^2 in red (•■•)• Since the vertex 25 (= 2n — fc) is saturated under Ai , we should ap- 
ply the injection In ■ In this case, X = {8, 9, 15, 16, 17, 21, 22} , cc = 20 and F = {1, 4, 5, 10, 13, 16, 21, 22, 23} . 
Since 20 is the sixth smallest element in X U {x} , we should choose 16 for y , which is the sixth smallest 
element in Y . Now cut the blue edge {20, 25} and join 16 and 24 with a red edge. From this, we can ex- 
tract two 8-matchings in K24, : ^i = { {1, 2}, {3, 4}, {5, 10}, {6, 11}, {7, 12}, {13, 14}, {18, 19}, {23, 24} } and 
B2 = { {2, 3}, {6, 7}, {8, 9}, {11, 12}, {14, 15}, {17, 18}, {19, 20}, {16, 24} } . (See FigureH) 

Remark 2. Theorem l4 . 1 1 and l4 . 21 iniplv the matching numbers {m{n, fc)} are log-concave for both coordinates: 

m{n, fc — 1) • m{n, fc + 1) < m{n, fc) for all fc > 1 , 

m{n — 1, fc) • m{n + 1, fc) < m{n, fc) for all n > 1 . 
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Figure 4: The injection In ■ 
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